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The interlayer magnetoresistance of layered metals in a tilted magnetic held is calculated for 
two distinct models for the interlayer transport. The first model involves coherent interlayer trans- 
port and makes use of results of semi-classical or Bloch-Boltzmann transport theory. The second 
model involves weakly incoherent interlayer transport where the electron is scattered many times 
within a layer before tunneling into the next layer. The results are relevant to the interpretation of 
experiments on angular-dependent magnetoresistance oscillations (AMRO) in quasi-one- and quasi- 
two-dimensional organic metals. We find that the dependence of the magnetoresistance on the 
direction of the magnetic field is identical for both models except when the field is almost parallel 
to the layers. An important implication of this result is that a three-dimensional Fermi surface 
is not necessary for the observation of the Yamaji and Banner oscillations seen in quasi-two- and 
quasi-one-dimensional metals, respectively. A universal expression is given for the dependence of 
the resistance at AMRO maxima and minima on the magnetic field and scattering time (and thus 
the temperature). We point out three distinctive features of coherent interlayer transport: (i) a beat 
frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angular- 
dependent magnetoresistance when the field is sufficiently large and parallel to the layers, and (iii) 
a crossover from a linear to a quadratic field dependence for the magnetoresistance when the field 
is parallel to the layers. Properties (i) and (ii) are compared with published experimental data for 
a range of quasi-two-dimensional organic metals. 
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I. INTRODUCTION 

One of the most fundamental concepts in electronic transport theory for elemental metals and semiconductors is 
that electronic transport involvesrthe coherent motion of electrons in band or Bloch states associated with well-defined 
wave vectors and group velocitieslil. An important issue is whether this roncept is applicable io interlayer transport in 
strongly correlated electron systems such as iigh- Tc superconductorsmfl, organic conductorsQ, and layered manganite 
compounds with colossal magnetoresistancela. If the interlayer transport is incoherent the motion between layers is 
diffusive and it is not possible to define band states extending over many layers and a Fermi velocity perpendicular to 
the layers. In that case a three-dimensional Fermi surface cannot be defined and Bloch-Boltzmann transport theory 
cannot describe the interlayer transport. I-I 

Extensive experimental studies have beenjirade of the angular-dependent magnetoresistance oscillations (AMRO)lI 
which occur in layered organic conductorsQ when the direction of the magnetic field is varied. The theoretical 
interpretation of these oscillations often involves a three-dimensional Fermi surface and their observation is sometimes 
interpreted as evidence for the existence of a thrce-dimcnsiopal Fermi surface. In quasi-|aiie-dimensional metals these 
effects are are known as DannerQ, Lebed (or magic angle )E3 E3, and third angular effectdiS, depending on whether the 
magnetic field is rotated in the a — c, b — c, or a — b plane, respectively. (The most- and least-conducting jiirections are 
the a and c axes, respectively). In quasi- two-dimensional systems, the effects observe include the Yamajil^-ascillations 
and the anomalous AMRO in the low-temperature phase of a-(BEDT-TTF)2MHg(SCN)4[M=K,Rb,Tlpla. 

The explaination of the Lebed effect is controversial and a number of different theories have been proposed.EHi^ 
It is not clear that coherent transport models can explain-lhe angjpjdependent magnetoresistaiice in the quasi-one- 
dimensional (TMTSF)2PF6 at pressures of about 10 kbarEiyiBOEa or the anomalous AMRO.N Consequently, we 
focus on the Banner and Yamaji oscillations here because their explaination in terms of a three-dimensional Fermi 
surface has generally been accepted. The resistance perpendicular to the layers is a maximum when the field direction 
is such that the electxoa velocity (perpendicular to the layers) averaged over its trajectories on the three-dimensional 
Fermi surface is zercoo. 

Several different models for incoherent interlayer transport have been considered previously. We shall distinguish 
between what we shall refer to as weakly and strongly incoherent interlayer transport. The former occurs when there 
is direct transfer of the electron from one layer to another and the intralayer momentum is conserved in the process. 
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Consequently, interference between the wavefunctions on adjacent layers is possible. However, the transport can be 
incoherent in the sense that tunneling events are uncorrelated because the electron is scattered many times W|itiiii3| 
the layer between tunneling events. This model has been used to describe interlayer transport in the cupratesEirE^I 
and organics.Eil In contrast, strongly incoherent transport occurs if the the intralayer momentum is not conserved 
by tunneling and there is no interference between the wavefunctions on adjacent layers. This can occuiJpecause the 
tunneling is associated with elastic scatteringp3, inelastic prcwjesses such as coupling to a bath of phononscJ or because 
of non-Fermi liquid effects such as spin-charge separation.Eli23 For both weakly and strongly incoherent transport, 
the interlayer conductivity of a bulk sample is determined by the tunneling rate between two adjacent layers. The 
interlayer resistance is then equal to the number of layers in the sample times the resistance between two layers. 

In this paper we present the details of calculations of the angular dependence of the interlayer magnetoresistance 
for both coherent and weakly incoherent interlayer transport when there is a Fermi liquid within each layer. Our main 
result is that coherent interlayer transport is not necessary to explain the Yamaji and Danner oscillations. Hence, 
their observation is not evidence for the existence of a three-dimensional Fermi surface. In contrast, we point out 
three properties of the interlayer magnetoresistance which do only occur if the interlayer transport is coherent: (i) 
a beat frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angle-dependent 
magnetoresistance when the field is parallel to the layers for sufficiently high fields, and (iii) a crossover from linear 
to quadratic field dependeace when the field is parallel to the layers. A brief report of some of the results presented 
here appeared previously.E3 

In the next section we present our main result, Eq. (|l]), an analytical expression for the interlayer conductivity in 
the presence of a magnetic field which is tilted at an angle 6 relative to the normal to the layers. This result is valid 
for incoherent transport for all field directions and for coherent transport provided the field is not almost parallel to 
the layers. We then use this expression to explain the basic features of the Danner and Yamaji oscillations. Simple 
expressions are described for the dependence of the interlayer resistance on the magnitude of the magnetic field and 
the scattering rate when the angle 9 is at an AMRO maxima or minima. In Section HI we derive Eq. (|^) for the 
case of coherent interlayer transport for both quasi-two- and quasi-one-dimensional systems. This involves evaluating 
Chambers' formula, a result of Bloch-Boltzmann transport theory. In Section ^ we derive Eq. for weakly 
incoherent transport in both quasi-two- and quasi-one-dimensional systems. In Section ^ we consider unambiguous 
signatures of coherent interlayer transport and compare these signatures to published experimental data for a range 
of quasi-two-dimensional metals. 



II. ANGULAR DEPENDENT MAGNETORESISTANCE OSCILLATIONS 

We assume that each layer of the metal is a Fermi liquid whose elementary excitations are fermions with wave 
vector {kx,ky) and with a dispersion relation of the form elkxyky). We consider the simplest posssible dispersion 
relations for quasi-one- and quasi-two-dimensional systems. (For a summary of our notation see Table I in Ref. p7| ). 
The interlayer conductivity in the absence of a magnetic field is denoted cr^z- We will show in this paper that in a 
tilted magnetic field the interlayer conductivity, for both coherent interlayer transport (except for fields very close to 
the layers) and weakly incoherent interlayer transport (for all field directions), is 



Jo(7tan6')^ 



(7 tan 6*) 2 
+ (i'WoTcos6')^ 



(1) 



where J^{x) is the v-th order Bessel function, loq is the oscillatioa, frequency associated with the magnetic field, and 7 
is a constant that depends on the geometry of the Fermi surface.EII The scattering time r is assumed to be independent 
of the momentum of the electron but can vary with temperature. 

If the field is sufhciently large and the temperature sufficiently low that loqt 3> 1 then the first term in (|l|) is 
dominant. However, if 'ytanO equals a zero of the zero-th order Bessel function then at that angle azz will be a 
minimum and the interlayer resistivity will be a maximum. If "/tanO 3> 1, then the zeroes occur at angles 0„ given 

by 



7tan0„ = TT(n " 

4:' 



(" = 1,2,3, 



(2) 



This condition was first derived for the quasi-two-dimensional case by YamajiQ and for the quasi-one-dimensional 
case by Danner, Kang, and ChakinEI. Determination of these angles experimentally provides a value for 7 and 
thus information about the intralayer Fermi surface. The values of the Fermi surface area of quasi-two-dimensional 
systems determined from AMRO are in good agreement with the Fermi surface areas determined from the frequency 



2 



of magneto-oscillationsB. Furthermore, AMRO can be used to map out the actual shape of the Fermi surface within 
the layer (see for example References ^,^8 29 ) . |— , 

The angular dependence of the interlayer resistivity p^z — ^/^zz^, given by Eq. (|^), for parameter values relevant 
for typical quasi-two-dimensional systems is shown in Fig. 0. Fig. |^ shows the angular dependence of pzz for 
parameter values relevant to (TMTSF)2C104. The results are similar to the experimental results in Ref. except 
near 90 degrees. Both figures are very similar to the results of numerical integration of Chambers' formula for 
coherent tcansport (Eq. (p7|)) except near 90 degrees. For coherent transport there is a small peak in Pzz{0) at = 90 
degrees. Q'EJ This is due to the existence of closed orbits on the Fermi surface when the field lies close to the plane of 
the layers. For incoherent tranport these orbits do not exist and so the associated magnetoresistance is not present. 
Since Eq. (|^) is also valid for incoherent interlayer transport the Banner and Yamaji oscillations can be explained 
equally well in terms of weakly incoherent transport. 



A. Asymptotic Form 

We want to find an expression for o■^^(^^) as 6* ^ §. Using the asymptotic formil 



which is valid for z >> , we can simplify Eq. (|l]) for the conductivity. Re- writing it as 



CTzz = 



Jo(m)' + 2 5]- 



^^1 1 + (^^)' 



(3) 



(4) 



where p, — 7tan0 and x = cjqt cos 9. We can substitute Eq. for J„(/i) provided that a; > 1 so that the sum 
in Eq. (^) converges rapidly. Separating the sum into the sum of even and odd terms, and using the fact that 
2 cos^ (a^ ~ ~ z) = 1 + sm{2p) for even n and 1 — sin(2/x) for odd n, gives 



TVp 



(1 + sin2;.) ( - + 5:^^^ ) + (1 - sin2,) ( ^.^^^ 



(2r 



(5) 



These series can be evaluated using the residue theoremc3 to give 



1 

px 



, , 7r\ sin(2/i) 
coth - ^ ' ^ 



sinh(f; 



(6) 



For ujqt > 1, Eq. (^) actually turns out to be a good approximation for 7 tan 6' > 1 (see Fig. 1). It will now be used 
to analyse the field and temperature dependence of the AMRO maxima and minima. 



B. Field and temperature dependence of the resistivity at critical angles 



Resistivity maxima: For the resistivity pzz to be a maximum 9 ~ 9^^^^ where 



7 tan 6*;; 



n - 1 I TT 



From this we can simplify sin(2/i) in Eq. (^) giving 

sin(2^) = sin(27tan6') = — cos(2n7r) = — 1 
for all n. The resistivity is then written as 

Pzz{OmaxiB) _ UJqT COS 



Pzz{B ^ {n - \) TT tanhf.^ 



COS 0„ 



(7) 



(8) 



(9) 
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This expression is plotted in Fig. |[ Now if the field is sufficiently high and the temperature sufficiently low that 
ujoT cos 6^g_^ >> 1, then the resistivity becomes 



-^(c^oT)^sin(2CaJ . (10) 



P,,{B = 0) TT 

Hence, at a fixed field the resistivity at the AMRO maxima will have the same temperature dependence as the 
scattering time which is inversely proportional to the zero-field resistivity. 

Resistivity minima: Similar arguments will show that for the interlayer resistivity to be a minimum 6 = 0'^^^ where 

7tanCm=U+7V (H) 



4 



and the resistivity then becomes 



p,4i? = 0)(n+i)7r ™" \2ujoTCOs 

This is plotted in Fig. ^. When luqt cos O^nin >> 1 then 

Pzz {^min ) _ 



^a;orcosC.ntanh . (12) 



+ - (13) 
p,,(B = 0) 2 V 4; ^ ^ 

which is independent of the field and scattering rate. Thus, at the AMRO minima the resistance will have the same 
temperature dependence as the zero-field resistance. 

Field in the layers: Now as — > ^, — > 71^0''", a; — > 0, and taking these limits in (^ gives 

The resistivity is linear in field at moderately high fields. However, caution is in order because in deriving (^) above 
we required that x > 1. That this is more restrictive than need be is suggested by the fact that Fig. 1 shows that 
remaias valid near 90 degrees. ladeed, Eq. ( p^ is valid: it agrees with the calculations of other authors for both 
coherentc3 and incoherent transportcJ in a quasi-two-dimensional system, provided the field is not to large (see Sections 
IV B and V C). Such a linear interlayer magnetoresistance has been observed in Sr2Ru04 (Ref. ^6|) and (TMTSF)2C10* 



(Ref. 37). However, for coherent transport the dependence on field becomes quadratic for sufficiently high fieldso 
(see Section 0). Note that ( |l^ ) is actually independent of the scattering time t. This means that for moderate fields 
the interlayer resistivity will only depend weakly on temperature. This was observed in (TMTSF)2C104: when the 
temperature was increased from 0.9 K to 8 K the zero-field resistivity (which is proportional to the scattering rate) 
increased by a factor of about six but the resistance at 12 Tesla increased by less than ten per cent. 

III. COHERENT INTERLAYER TRANSPORT 

If the interlayer transport is coherent one can define a wave vector perpendicular to the layers and a three- 
dimensional dispersion relation esD^k) of the form 

esD {k) ^ e{kx,ky) - 2tc cos{k^c) , (15) 

where tc is the interlayer hopping integral, c is the layer separation, and e{kx, ky) is the intra-layer dispersion relation, 
simple examples of which are given below. The electronic group velocity perpendicular to the layers is 

lde{k) 2ct, 

.. = ^^ = — sm(cM. (16) 

We calculate the interlayer conductivity by solving the Boltzmann equation in the relaxation time approximation 
which leads to Chambers' formulafil 



df{E) \ 



= Mk)vAk) [--^ ) d'k , (17) 
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where /(e) is the Fermi function and r is the scattering time which is assumed to be the same at all points on the 
Fermi surface. Vz{k) is the electron velocity perpendicular to the layers and is averaged over its trajectories on the 
Fermi surface 

Vz{k) = ^J exp (J-^ v4k{t))dt . (18) 

where k{0) = k. The time dependence of the wave vector k{t) is found by integrating the semi-classical equation of 
motion 



dk e ^ -:■ , „, 

37 = -72 Vfe^sD X B . (19) 
at Ti 



Now if the temperature is sufficiently low that T « Ep then ^ in Eq. (17) can be replaced by a delta function at 
the Fermi energy and Eq. (|l7|) becomes 



47^3 



v,{k)v,{k)5{EF - eMk))d''k . (20) 



A. Quasi-two-dimensional case 

Here we consider a quasi-two-dimensional system with the energy dispersion relation 

^2 

eMk) = ^ (fc' + ^y) - cos (fc,c) , (21) 

where m* is the effective mass of the electron. We assume the interlayer hopping is sufficiently small that tc <C Ep- 
The Fermi surface is then a warped cylinder (see Ref. ( p7|)). Substituting the energy dispersion relation from Eq. ( [2l| ) 
we obtain the components of the group velocity 

1 -± f hkx hky 2ctcSui{ckz) 



m = ^r^/^, — 7^^-- . (22) 

n \ m m a / 

In order to calculate the time dependence of kz we must integrate Eq. ([l^ ) which can be written in the form 

^ — {—kyB cos9,kxB cos9,kyB sinO) . (23) 
dt m* 



Terms of order tc have been neglected once we assume tcta.n9 <C , where kp is the Fermi wave vector, defined 

by Ep = . Differentating the x and y components of Eq. (DSh with respect to time we obtain a second order 
differential equation whose solution gives ky[t) — kp cos{uJct) and kj.{t) — kp sai{LOct), and 

eB cos 6 , , 

— 24 

m* 

is the cyclotron frequency. Substitution of this into the z component of Eq. ( p3| ) and integrating gives 

A:^(i) = fc^(O) + A:i?tan6'sin(wci) . (25) 

In order to calculate the z component of the group velocity we substitute the expression for fc^ (t) into the z component 
of Eq. dH) giving 

Vz(kz{^)A) = Csin(Msin0 + fc^(O)c) , (26) 

where 

/i = ckp tanS , (27) 
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(f> is the angle around the orbit, and C = 2ctclh. Integrating the velocity in Eq. (^6|) over a period gives us the average 
velocity which can be written as 

< Vz >oc / sin(/x sin Jo(7tan0) (28) 

^0 

and in the absence of scattering this average velocity is equal to zero when jtanO equals a zero of Jq. These particular 
values of 9 corrrespond to the peaks in the resistivity. 

We can write Eq. ( p^ in a slightly simplified form in order to highlight the fact that the integral is a surface integral. 
If the warping of the Fermi surface is small we can parametrise the surface using kz and 4> where kx = kp cos (p and 
ky = kp sin (j) giving 



dS Vz{k)vz{k) 



e rm ' ' 



FS 



47r3r 



dk. 



tt/c 



d(f>Vzik)vz{k) 



(29) 



here Vz{k) is defined in Eq. ( |l8|) and the pre- factor m* /h^ arises from the delta function . In terms of the parametrised 
surface we have 



Vz{k)^ / exp((/)7™c) i'z(fcz(0), - (/)') 

- TiOr 



(30) 



where we used the fact that = ujct and Vz{k{t)) = Vz{kz{Q),<j> — (f>'). For closed electron orbits the electron group 
velocities are periodic functions of (j) and ((>' . Thus the range of integration of (j)' can be cut up into segments each 
having length 27r. The conductivity is then 



4^ 



7r/c ^ 

dkziO)- — — 

-tt/c (1 - exp(-27r/TWc) 



d(l)Vz{kziO),(f>) X 











Vzikz{0),(j)-cb')exp{(b'/Tu;c) 



l-2n 

We use trigonometric identities to expand Eq. (^) and substitute the Bessel generating functionsH to obtain 



(31) 



Vz{kz{0),q^~cb')^C sin(fc,(0)c) 

C cos{kz{0)c) 
If we substitute this into Eq. (pO|) we obtain 



Jo(m) + 2 J2 JM cos((2A:)(0 - ^')) 



k=l 









(^sm{kziO)c) 


J-oo TUJc 



2 J2k+i{fi) sm{{2k + 1)(0 - (/)')) 
. fe=o 



Jo(m) + 2 J2 cos((2fc)(0 - 0')) 



(32) 



k=l 



cos{kz{Q)c) 



2 J2k+i{fi) sin((2fc + - cj)')) 



L k=0 



exp((/) /tlOc) 



(33) 



Substituting equations for Vz{kz{0), (f) — (j)') and Vz{kz{0),4>) into Eq. ( ^l| ) we note, terms that survive is when k = l, 



since integrals such as /^^^ c?(?!)cos(2fc0) cos(2/(/)) = ttSm, where 5ki is the Kronecker delta, thus giving 



dkz{0) [2TrTUJcSm{kz{0)cfM^^f 

J2k{piY 



^.r ^-^(^^W^)'E-(2fc)^ + (l/.,.)2 



cos(fc,(0)c)2^ 



J2k+l{l^y 



k=0 



(2fc + 1)2 + (l/weT)^ 



(34) 



Performing the integral over ^^(0) yields the feial expression for the conductivity which is of the form of Eq. (1) 
This result was previously given by Yagi et alsB. 
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B. Quasi-one-dimensional case 



For the quasi-one dimensional case we begin with the dispersion relationH 

^3D{k) — hvF{\kx\ — kp) — 2tb cos{kyb) — 2tc cos{kzc) (35) 

where vp is the Fermi velocity and tb is the intrachain hopping within the layers. The Fermi surface consists of two 
sheets at k^ = ikp. By proceeding as for the quasi- two-dimensional case the rate of change of the wave vector (and 
defining the magnetic field B ^ {B sin 9, 0, B cos 9)) is given by 



^ / —2bB cosdett,sm(bky ) 



eBh cos OvF \ b ] , (36) 

26_Bsin6'eif, sin(6fcj,) / 



where we neglect terms involving tc by assuming that vp ^ Vz tanO, i.e., the warping of the Fermi surface is small 
and the magnetic field is not too close to the layers. Integrating the second equation, gives 



ky{t)^'^t + ky{0) , (37) 



eBbcos9vp 

ujB = (38) 



where 



is the speed at which the wave vector traverses the Fermi surface. To obtain kz{t) we substitute Eq. (|3^) into (|36|c) 
and integrate to obtain 

. / s , / ^ 2ebthB sin 9 , , , 

kz{t) =fc^(0) \ cos{ujBt + bky{Q)) . (39) 

Substitution into the z-component of the velocity yields 

, 2ct / 
?;,(fc,(O),0-0) = -^sin(cA:,(O)-7tan0cos(0-0)) , (40) 

where = —WBt, (j) = bky{0) and 7 = This is similar in form to the z component of the velocity for the quasi-two 
dimensional case (compare with Eq. (|^)). The interlayer conductivity can then be written as 



azz^-^jT— I dkz{0) dcl>Vz{kz{Q),<l>) 



— n -T^i ^^^.(fc.(O),0-0) • (41) 

(1 - exp(-27r/rwB)) 

The integral from —2tt to over is obtained by noting that the electron group velocity is a periodic function of 
(f) , we can cut the range of (p into segments having length 27r. The factor (1 — e^^'^^'^"^) is a consequence of this. 
Proceeding as for the quasi-two dimensional case leads to the result of the form Eq. (^. As far as we are aware this 
expression has not been derived previously for quasi-one-dimensional systems. 



IV. WEAKLY INCOHERENT INTERLAYER TRANSPORT 



Suppose that the coupling between the layers is sufficiently weak that the time it takes an electron to hop between 
the layers (approximately fi/tc) is much longer than the scattering time. This means that the intralayer scattering 
rate is much larger than the interlayer hopping integral 

- > (42) 
r 

and the mean-free path perpendicular to the layers is much smaller than the interlayer spacing. If this condition 
holds then the interlayer transport will be incoherent in the sense that successive interlayer tunneling events are 
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uncorrelatecS. Previous estimates of tc and r in various layered organic metalJ^'BEl suggest these quantities may be. 
comparable at low temperatures. Furthermore, the scattering rate usually increases quadratically wi|tli temperatureuHl 
and at temperatures of the order of tens of Kelvin this condition will almost certainly be satisfiedE^ The interlayer 
conductivity is then proportional to the tunneling rate between just two adjacent layers. If we assume that the 
intralayer momentum is conserved the tunneling rate can be calculated using standard formalisms for tunneling in 
metal-insulator-metal junctions.c3 Modelling the interlayer transport in this way is reasonable because many organic 
conductors consist of conducting layers separated by insulating layers of anions that are several A thick. Furthermore, 
intrinsic Josephson type effects have been observed in the superconducting state of k-(BEDT-TTF)2Cu(NCS)2 (Ref. 

. . . 

We consider the simplest possible model for the tunneling between layers, direct transfer described by the Hamil- 
tonian 

Hi2^-tcj d?r (ci(r)tc2(r)+C2(r)tci(r)) , (43) 

where ci(r)^ creates an electron in layer 1 at r. Note that the interlayer transport is coherent in the sense that during 
the tunneling process the phase information in the electron's wave function is not completely lost. However, it is 
incoherent in the sense that due to the large intralayer scattering rate the interlayer transport cannot be described by 
Bloch states extending over many layers. If we consider a sequence of tunneling events they are uncorrelated because 
after a tunneling event an electron is scattered many times before it tunnels to the next layer. 

The interlayer current / associated with iJ]|a-.and produced by a voltage V can be calculated using the formalism 
developed for metal-insulator-metal junctions.uS The result for the interlayer conductivity is 



dl 



L^Ly dV 



y^=nLi-J " n " ' J ^^^(^''^ '^)^^(-''^ '^)^ ' (44) 



where Ai and A2 are the spectral functions for layers 1 and 2 and and Ly are the dimensions of the layers. It 
will be seen below that in the presence of a tilted magnetic field Ai and A2 are not identical. The zero-field limit 
(for which Ai =_4glx>f this expression has been used in treatments of incoherent interlayer transport in the cuprate 



superconductorsc3E3l23. 



If we assume T << Ep , then can be replaced with a delta function to give 



OE 



d'r / d'r Ai{r,r\EF)A2ir\r,EF) . (45) 



This can be re-written by noting that the spectral fimction can be written as Ai^2{f' Ep) = —^[6*^2(^1 ^'j ^f") ~ 
^r,2(^' ' Ep)] leading to 

<7..{r,f, Ep) = J d'r J d\' [G+(f, r\Ep)G2 {P ,r, Ep) + G^ir\r, Ep)G+{f, P,Ep)] . (46) 

In the Landau gauge, the vector potential A, for the magnetic field B = {Bx,0, Bz) = {B s'm9, 0, B cos6), is 

A=iO,xBz-zB,,0) , (47) 

where B and A are related by _B = V x A. The vector potential in the two layers (see Ref. ( |2^)) are not equal and 
differ by a gauge transformation A2 = Ai + VA where 

VA = Ai - i*2 = (0, ~cB sin 9, 0) . (48) 

The Green's functions in the two layers are not equal. This reflects the fact that even though the magnetic field is 
invariant under a gauge transformation the Green's function is not. The Green's function for layer 1 is thus multiplied 

by a phase factor exp ||-e[A(r) — A(r')]| giving 

G+(f, /) =exp|jeA(f)| G^(r,/) exp|^eA(/)| . (49) 
Making use of this relationship we have 
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/ \G+{r,Q,EF)\\os{^sm0y) , (50) 



where we have used the fact that ICf (r, 0, Ep)\'^ is translation invariant. 



Note that Eq. (50) is a very general expression which holds provided that intralayer momentum is conserved and 



the interactions between, the layers can be neglected. It is valid in the presence of interactions within the layers and 



for a non-Fermi liquid.Ej Second, this expression shows that for weakly incoherent interlayer transport the intcrlayer 
conductivity is completely determined by the_,one-electron Green's function whereas the intralayer conductivity is 
determined by two-electron Green's functionsEJ. 



It is the averaging of the phase factor over the spatial integral in Eq. (50) that gives rise to the AMRO effect. 
The length scale associated with the magnetic field for the quasi-two-dimensional system is the cyclotron length R 
which at the Fermi energy is i? = Tikp / {eB cos 0). For the quasi-|ape-dimensional case the length scale associated with 
oscillations perpendicular to the chains is i? = 2%/ [evpB cosO)^. At this length scale the phase difference between 
the wave function of adjacent layers is A(i?) = eB smOcR = ^ianO. Naively, we might expect maximum resistivity 
when this phase difference is an odd multiple of tt, leading to a condition different from Eq. (^). However, one must 
take into account averaging of the electron position over the perpendicular direction. 

We now proceed to evaluate (|5^) for the simplest possible situation, where there is a Fermi liquid within each layer 
and the magnetic field is small enough that we can take the semi-classical limit of the Greens functions. 



A. Quasi-two-dimensional case 
The Green's function for layer 1, in the absence of scattering can be writterS 



G(f,r',<)i- 



ct/2 



2T:iht sin(wct/2) 



exp 



-L 



and 



L = 



^ cot(^) + (x + x )(y -y) 



(51) 



(52) 



where lDc = io^cosd = eBcoaO/m* is the cyclotmn frequency. In order to calculate the conductivity we follow the 
same approach that Hackenbroich and von OpperO used to study Shubnikov - de Haas oscillations in two-dimensional 
electron systems. In the presence of scattering the energy-dependent Green's function is 



(53) 



G+{f,r' ,E) = — I dtexp{ -{E + iT)t] G{r,r' ,t) 



where F = ^ is the scattering rate. The retarded Green's functon is obtained using G~{r' ,r,E) = [G~^ {r, r' , E)]* . 
We perform the integral in Eq. ( p3| ) by the stationary phase method which is valid in the semi-classical limit (?i 0). 
The stationary phase condition gives 



E = 



m*ujl 



sa\{uict/2) 



(54) 



This shows that if the cycloron radius is Rc then G^{f,r' ,E) vanishes for \r — r'\ > 2Rc, while, jQr |r — r' | < 2Rc 
there exists two different cyclotron orbits and one finds an infinite set of stationary times given bjc3 



2nn 



(55) 



where n determines the number of revolutions the electron makes to get from r to r' and q = S or L denoting the 
two different paths it can take (see Fig. ||) . The times to traverse these paths are calculated using the the stationary 
phase condition and E — m* R^uj'^/2, to give 
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Putting all this together and performing the integrals we obtain 

G+{f,P,E) 



m 
2ih 



EE 

n=0 q=S,L 



exp 



T 



n.q 

2^ 



I ITT 
X exp [ —On.q '7rVn,q 



where 



^n,q — ^T^n.q 



2 cot( — Y'^) + ix + x ){y -y) 



(57) 



(58) 



and 77 is the Maslov or Morse index (the number of conjugate points along the orbit (Ref. |46[ p. 223)), r/n^s = 2n and 
Vn,L = 2ri + 1. Equation (pQ) can be written as 



d^\r — r' \ 



G^{r,r,EF) ^ cos(^Bsin6'| 



r — r sm < 



(59) 



where cj) is the angle between the vector \f — r'\ and the x-axis. Substituting the semiclassical expressions for the 
Green's functions into the above equation and changing the integrals over r and r to polar coordinates, one obtains 
a double sum (denoted by the subscripts (1,2) of n and q) over the classical trajectories 



ni,Tl2=0 qi,q2 = S,L 



[Tni,qi + Tn2,q2] 



2r 



^sin(wcr„i^gj sin(wcT„2,g2) 
cos 



2, 92 



■ smy r — r smi 



(60) 



We make the simplification that qi = q2, for when qi ^ q2 the integrand oscillates, that is as ?i ^ the oscillations 
cancel each other and therefore do not contribute. This gives 



etc\ f LOcC 



h\ J \2Ef 



/ dcj) I \r — r' \d\f — r' \ > exp (ni 

Jo Jo ...tli'-n V ^cT 



"2 



7ii,n2— 



exp(--ts/r) exp(-tL/T) 



sm{ujcts) sm{uJctL) 



cos 27r 



Ef 1 
TitJc 2 



(n2 - "^i) I cos 



f ecB 

SI 



smW r — r sm( 



(61) 



Terms with ni ^ n2 correspond to the Shubnikov-deHaas oscillations. We neglect these by setting ni = n2 since they 
will be smaller that the leading order terms by a factor of order exp( — —) and thus we have 



2ix ju. 



exp 



7r/ci 



T 



cos r\ sm 



exp 



is 

T 



COS 77 sm 



(62) 



where A — and 77 = ^j^B svaORcSaKf) = 2cfcF tan 6' sin — 27tan6'sin(/). Combining the integrations over ts 

and II and performing the summation over n, one obtains 



10 





p27r/ujc 




/ exp(— i/T)cos 


Jo 


Jo 



k=l 



To evaluate the integral over t we raake use of the identityS 

oo 

cos(rysin/3) = Jo(r/) + 2^ J2fe(?/) cos(2fc/?) 
The conductivity then simplifies to 

CTzz =At Jo(77) + 2 V- 

•^0 L fc=i ^ 



This integral is of the form j^^^ J2kiz aiii(j))d(j), where z = 2^tan0, which can be evaluated using the relatior 



(63) 



(64) 



(65) 



/ J2A: (z sin 0)ci0 = 27rJfc (z/2) 
Jo 

conductivity wl 

interlayer tunneling of a quasi-two-dimensional system in the absence of scattering! 



(66) 



We then obtain an expression for the conductivity which is of the form Eq. (|^). Pceyiously, Yoshioka calculated the 

rinK.Ea 



B. Quasi-one-dimensional case 

The Hamiltonian within a layer in a magnetic field is 

d f d 

H — aifiVF-^ 2ihC0S b -— exBcosO 

ox \ \i oy 



(67) 



where a = ±1 denotes which sheet of the Fermi surface the electron is on. The wave function within a layer is given 
byEJ 



i^k^,ky,a{x, N, t) = exp \ *[-^ + + bkyN - aXsui{kyb - qx)] 



(68) 



where N denotes the number of the chain, x is the distance along the chain and the dispersion relation 

eaikx, ky) = ahkxVF (69) 

and 

ebB cos ujb 



q = 



h 



VF 



where lub is the oscillation frequency given by (pq) and 



A 



2th 



ebvFB cost 



(70) 



(71) 



The transverse motion of the electrons due to the field is approximately XbsB The one-electron advanced Green's 
function in the absence of scattering is 



G+ix,x,N,N',t,0)^ rk^,k,A^',N,t)i^k^,k,Ax,N,0) 

kj;,ky,CX 

for t > 0. Taking the Fourier transform (with respect to time) of this and including a scattering rate ^ — jip 

1 



(72) 



G+(a;,a;',iV, N , E) 



^ , dteyii>{-{E + iT)t'> G+ {x, x , N ,t,0) 
in ./o la 



(73) 
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After performing the integral in t and 



G+ = 



^ E - N') + aXL] } exp (fc^ (e 



Tivp ^ — ' 01 

ky ,Q 



2t 



where 



L = sin{kyb — qx ) — sin{kyb — qx) . 

This is similar to the quasi-classical Green's function given by Gorkov and Lebed@ 
The conductivity (Eel) then becomes 



(74) 



(75) 



X expjiaAS'i} exp 



dx / dx 



exp <^ i{N - ) kyib - ky2b 



ebcB sin 









> + c.c 


y VFT 





where 

Si = Li — L2 = sm{kyib — qx) — svnikyib — qx) — sin(fcj,2& ^ ix) + sin(fcj,2& — qx) . 

If we now let M± ^ N ± N then the sum over M_ gives a deha function 2TrS{kyib, ky^b - 
with kyb the conductivity then simplifies to 



— dx dx 



exp{iQ;AS'2} exp 



f \x-x'\ 






> +C.C 







(76) 
(77) 

Replacing kyib 
(78) 



where 



S'2 = sin(fcj,6 ~ qx ) — s\TL{kyb — qx) — sin kyb 



ebcB sin 6 



qx 



sin I fcj^6 — 



ebcB sin 



(79) 



Now shift kyb to kyb + qx and make the substitution A 



ebcB sin 6 



-j^ , thus 52 becomes 

S'2 — sm{kyb) — sin(fcj,6 — q{x — x )) ~ sm{kyb — A) + sin(A:y6 — A — q{x ~ x )) . (80) 
We now let x± = ^ and a;_ — vpt and perform the integral over x^. This simplifies the conductivity giving 



/•oo 

/ exp { I a A exp 



2t 

T 



+ C.C , 



(81) 



where ^3 is given by 

5*3 = sin(fcj,6) — sin(fcj^6 — ^ust) — sm{kyb — A) + sin(fey& — A — 2u)Bt) (82) 
and ljJb ~ qvp. This can be separated into two parts and simplified using the appropriate trigonometric identities 
83 = n- f3 



2 cos ( kyb — — I sin ( — | — 2 cos ( kyb — - — 2cjRt I sin ( — 



A cos ( k„b — ^ I — A cos ( kyb 2wst 



(83) 



where we have taken A << 1. We can justify this by considering the dimensions of the unit cell and the magnetic 
fiux passing through the area of the cell. If the magnetic field i? ~ 10 T, and the area be ^ lO^^'* m^, then the flux 
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$ = Bbc will be small and thus A = $/<I>o << 1 where <J>o = ^/e is a flux quantum. Re-writing the conductivity we 
obtain 



^ f°° , r 2t] r''' dky , , , 

> / aiexp< > / — -2- exp {— laAp} exp {laA/x} + c.c 

„ Jo I T J J-Wb 27r 



Using the identityEj 



exp 



-{h-k 
2^ h' 



the exponentials in ^ and p become 

oo 

exp{— zaA/3} — ^ — Jn(aAA) exp 

n— — oo 

oo ^ 

exp {iaXfi} — J^/ (aAA) exp < m 



n — --00 



A 

kyb - — ~ 2ujBt 
2 



Substitution of these into Eq. (p4|) and performing the integral in t gives 



oo oo 



^iS^E E E J«MA)J„,(aAA) 

a n=-oo„'^_oo 



1 



1 + inijJBT 

This integral is zero unless n — —n , thus 



-7r/b 



dfcj^ exp < I I kyb — ;7 I + c.c 



bh nvp 



a n— — oo 



J„(a2:)^ 



1 + inuJBT 



+ c.c 



(84) 



(85) 



(86) 



(87) 



(88) 



where z = ^ = % g = 7tan0 and 7 is the same as for the coherent case. The summation over a is performed 
by noting that Jn{z)'^ = Jn{—z)'^ foi' all n. Finally we include the complex conjugate to obtain an expression which 
can be written as Eq. (m. 



C. The magnetic field parallel to the layers 



We consider here the field range over which the result (14) holds for incoherent interlayer transport. We define 
the magnetic field as i? = {Bx,0, 0) and the vector potential as A = (0, —zB^, 0). It is easiest to work with spectral 
functions in the momentum representation; the interlayer conductivity is then given by 

= ^^^l(fc,i?F)v42(^,i?F) , (89) 

k 

where Ai and A2 are the spectral functions for the two layers. Due to momentum conversation we can then write the 
spectral function for layer two in terms of layer one as 

Ai{k,EF) = A2(k^^A,EF) ^A{k,EF) . (90) 



where 



2T 

(£;f -e(fc))2 + r2 



13 



and r = Ti/2t and e{k) is the dispersion within the layer. 

We now specialise to the quasi-two-dimensional case. Substituting (|9l|) and (BOl) into the conductivity gives 



Now, introduce polar co-ordinates {k,(f)) so k1 
becomes 



kl 



+ r2 

k 



Ep — 



2m* 



k'^r ^" ( ku 



(92) 



r2 



V 



kcoscj), and define A = [ecB/fif' so that Eq. (p2 



1 ^2ir 

Jo 



r 



.2 



r2 



^fc2 _ (k^ + A+'^cBkcosi 

2m* -T 2m* V ?i 



r2 



(93) 



Suppose that the field is sufficiently large that F ^ heEkpc/m* (which corresponds to ojqT'^ ^ 1) then the first 
spectral function has a sharp peak near k — kp whereas near that peak the second term varies slowly. Hence, we set 
k = kp in the second term and then integrate over k to give 



e^tjcTm* 



1 



{A +fcBkp COS. 



(94) 



r2 



When A <^ ecBkp/fi, this integral will be dominated by the behaviour near the two zeros of coscj) so we can write 
the integral as 



2eHlcVm* 



[±,ecBkpcf\ +r2 



(95) 



Performing the integral gives (|T^) resulting in a magnetoresistance which is linear in field. 

When A ~ ecBkp/h that is ecB/hkp ~ 1 deviations from this linear in field behaviour will occur. If c ~ 10 A, 
ckp ~ 3, then B « 2000 T. Similar arguments apply to the quasi-one-dimensional case. It will be shown in the next 
section that for coherent inter layer transport the deviations from linear dependence can occur at much lower fields. 

V. DEFINITIVE TESTS FOR COHERENT INTERLAYER TRANSPORT 

We have shown that the Yamaji and Danner oscillations exist for both coherent and weakly incoherent interlayer 
transport and so cannot be used to establish that the Fermi surface is three-dimensional. We now consider three 
properties which are different for coherent and incoherent interlayer transport. 



A. Beats in magnetic oscillations 



For quasi-two-dimensional systems definitive evidence for the existence of a three-dimensional Fermi surface, is the 
observation of a beat frequency in de Haas-van Alphen and Shubnikov - de Haas oscillationsHD The frequency F of 
these oscillations is determined by extremal areas A of the Fermi surface, F = hA/ {2Tre) {Re f. 0). For the warped 
cylindrical Fermi surface (see Fig. 1 in Ref. ^) there are two extremal areas, corresponding to the "neck" and "belly" 
orbits. The small difference between the two areas leads to a beating of the corresponding frequencies Fi and F2. In 
a tilted magnetic field the frequency difference is 



Fi — F2 
Fi 



-—^Jo{kpcta,nt 
Ep 



(96) 



Table | lists several materials in which such beat frequencies have been observed. In_/3-(BEDT-TTF)2l3 and /3- 
(BEDT-TTF)2lBr2— the angular dependence of the beat frequency is consistent with (|9q ) and tc/Ep ~ 1/175 and 
1/280, respectivel}!^!. 

However, Table ^ indicates that in many other quasi-two-dimensional organic metals no beat frequency has been 
observed. This could be because the interlayer transport is incoherent or because the interlayer hopping tc is so small 
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that the beats cannot be resolved experimentally. Suppose that oscillations but no beats are seen in the field range 
from Bmin to Bmax- This means that cos(27r(_F'i — F2)/B) has no zeroes in this field range, implying that 

ti~ t2 < — — ^ . (,97j 

This together with Eq. ( p6| ) can be used to establish apjipper bound for td Ep. For K-(BEDT-TTF)2l3 the absence 
of beating has been used to establish tc/Ep < l/3000l3o. This implies a resistivity anisotropy Pxx/Pzz (tc/Ep)^ < 
10~^. However, the observecS anisotropy in the k-(BEDT-TTF)2X materials is about 10~^. This inconsistency 
suggests that the interlayer transport may be incoherent in K-(BEDT-TTF)2l3. However, it could be that the measured 
value of 10^"^ is too large because resistivity anisotropy is too large because the measurement of pxx involves some 
component of p^z due to an imhomogeneous current distribution or the current path being changed by sample defects. 



B. Peak in the angle-dependent magnetoresistance at 90 degrees 

Numerical, solutions of Chambers' fopnula (^) for coherent interlayer transport show that for both quasi-one- 
dimensionalEI and quasi-two-dimensionalEll materials, at sufficiently high fields, the angle-dependent magnetoresistance 
has a peak as the field direction approaches the layers (i|e., at = 90 degrees). This peak is absent for incoherent 
interlayer transport (see Figs. |l| and Hanasaki et al.O identified the peak as being due to closed orbits which 
occur when the field is parallel to the layers. These orbits are associated with the cyclotron frequency 



n^u,[^-^] =u,^[^] (98) 



and so will only be important when the field is sufficiently large that fJr > 1. 

Table || lists whether or not the peak has been observed for a range of quasi-two-dimensional metals. Note that 
the presence (absence) of the peak is not always consistent with the observed presence (absence) of beats. This can 
be because the two sets of measurements were done on different samples of different purity (and thus had different 
values of r) or because the field was not large enough to observe the peak. The presence of a peak at 90 degrees in 
the AMRO dataO for (TMTSF)2C104 suggests that it has coherent interlayer transport. 



C. Crossover from linear to quadratic field dependence for a magnetic field parallel to the layer 



Schofield, Wheatley and Coope considered the interlayer magnetoresistance for quasi-two-dimensional systems 
with coherent interlayer transport and a magnetic field parallel to the layer. Eq. (25) of Reference 34 gives an 
expression for the interlayer conductivity for all values of the magnetic field. They showed that when U,t <^ 1 the 
magnetoresistance increases linearly with field, as in Eq. (14). However, for fir ^ 1 the field dependence becomes 
quadratic and is given by 



1.96 



m_ 

,(0) (7Wot)2 
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1/2 



(99) 



The deviations from linear behavior will occur when fir > 1, i.e, ujqt > i (■f^j • ^'^^ typical organic samples this 
will happen in the field range of 10-100 tesla. In contrast for incoherent interlayer transport it was shown in Section 
IV C that the deviation from the linear field dependence would not occur until about 2000 T. We are unaware of 
any material in which a search for this linear to quadratic crossover has been made. This field dependence is to be 
contrasted to that at angles slightly different from 90 degrees, which will be given by Eq. (p^. The ratio of these two 
expressions provides a means to determine tc/Ep since 7 and loqt can be deduced from AMRO data. 



VI. CONCLUSIONS 



We have presented detailed calculations of the interlayer magnetoresistance of quasi-one- and quasi-two-dimensional 
Fermi liquids in a tilted magnetic field. Two distinct models were used for the interlayer transport. The first involved 
coherent interlayer transport and made use of results of semi-classical or Bloch-Boltzmann transport theory. The 
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second model involved weakly incoherent interlayer transport where the electron is scattered many times within a 
layer before coherently tunneling into the next layer. We found that the dependence of the interlayer magneoresistance 
on the direction of the magnetic field is identical for both models except when the field is almost parallel to the layers. 
An important implication of this result is that coherent transport is not necessary for the observation of the Yamaji 
and Banner oscillations. Hence, observation of one of these effects in a particular material cannot be interpreted as 
evidence that the material has a three-dimensional Fermi surface. Instead, we propose three unambiguous tests for 
coherent interlayer transport: (i) a beat frequency in the magnetic oscillations in quasi-two-dimensional systems, (ii) 
a peak in the angular-dependent magnetoresistance when the field is parallel to the layers, and (iii) a crossover from 
a linear to a quadratic field dependence for the interlayer magnetoresistance when the field is parallel to the layers. 
A survey of published experimental data on a wide range of quasi-two-dimensional organic metals suggests that some 
have properties (i) and (ii) and others do not. 

In future publications we will examine the frequency dependent interlayer conductivity and the Lebed and 
third angular effects in quasi-one-dimensional systems. iL-|piuch greater challenge is to explain the AMRO ob- 
served in (TMTSF)2PF6 at pressures of about 10 kbaiMlla and in the low-temperature phase of a-(BEDT- 
TTF)2MIIg(SCN)4[M=K,Rb,Tl]. The angular dependence of the latter is inverted compared to that of the Yamaji 
effect. In particular, the magnetoresistance is smallest when the field is in the layers, the opposite of what one expects 
based on the simple Lorentz force arguments relevant to semi-classical magnetoresistance. Understanding this may 
require knowledge of the effect of an orbital magnetic field on a strongly correlated electron system. Little is known 
about this problem except in the extreme quantum limit of a partially filled lowest Landau leveS which is far from 
the situation considered here where usually of the order of tens of Landau levels are filled. 
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FIG. 1. Angular-dependent magnetoresistance oscillations. The dependence of the interlayer resistance for a typical 
quasi-two-dimensional system on the direction of the magnetic field is shown for a range of magnetic fields. 9 is the angle 
between the field and the normal to the layers, r is the scattering time within the layers and ljq is the cyclotron fr equency 
when the field is perpendicular to the layers. The curves shown are plots of Eq. (1) which is valid for all 6 for incoherent 
interlayer transport and for all 9 except close to 90 degrees for coherent interlayer transport. Note that the location of the 
maxima and minima is independent of the field and the scattering time. The dashed curve is a plot of the asymptotic expression 
(0), which can be seen to be a very good approximation for 9 > 20 degrees. 
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FIG. 2. Dependence of the interlayer resistance of a quasi-one-dimensional system on the direction of the magnetic field. 
8 is the angle between the magnetic field and the least conducting direction, with the field in the same plane as the most 
conducting direction. The parameter which defines the anisotropy of the intralayer hopping 7 = |^ = 0.25. r is the intralayer 
scattering time and is the frequency at which the electrons oscillate between the chains when the field is perpendicular to 
the layers. Except very close to 90 degrees this figure is similar to the experimental data on (TMTSF)2C104 in Ref. 0. 
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FIG. 3. Universal dependence of the interlayer resistivity on the magnetic field and scattering time when the field is tilted 
at an angle corresponding to an AMRO minimum {9min) and an AMRO maximum (Omax)- For high fields the resistivity at 
the minima becomes independent of field and has the same temperature dependence as the zero-field resistivity. For high fields 
the resistivity at the maxima increases quadratically with field and has the same temperature dependence as the inverse of the 
zero-field resistivity. The curves are not plotted for small ujot cos 9 because the results derived in the text are not valid in that 
regime. 
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FIG. 4. Short and long semi-classical orbits joining two points within a layer of a quasi-two-dimensional system in a magnetic 
field perpendicular to the layers. Rc is the radius of the cyclotron orbit. 
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TABLE I. For a range of quasi-two-dimensional materials we list whether or not beats in magnetic oscillations and a peak 
in the angular dependent magnetoresistance at 90 degrees has been observed. For coherent interlayer transport both these 
features should be present provided a wide enough range of magnetic fields is explored. A question mark indicates that the 
measurement has not been made. 

Poak at 90 degrees 



a-(BEDT-TTF)2NH4Hg(SCN)4 

a-(BEDT-TTF)2KHg(SeCN)4 

a-(BEDT-TTF)2KHg(SCN)4 above 20 T 

a-(BEDT-TTF)2TlHg(SeCN)4 

a-Et2Me2N[Ni(dmit)2]2 

a-(BEDT-TSF)2KHg(SCN)4 above 6 kbar 

/3h-(BEDT-TTF)213 

/3-(BEDT-TTF)2lBr2 

K-(BEDT-TTF)2l3 

K-(BEDT-TTF)2Cu2(CN)3 at 7 kbar 

K-(BEDT-TTF)2Cu(SCN)2 

e-(BEDT-TTF)2l3 

Sr2Ru04 



nc 
nc 
nc 

yeJ 

7 

yei 
yei 
yei 
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